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ABSTRACT 

Using the BeppoSAX observatory, we have observed a nearby LINER/Seyfert 2 galaxy, NGC 
3079, which is known as an outflow galaxy and a bright H2 0-niaser source. Using the PDS 
detector, we have revealed that the NGC 3079 nucleus suffers from a Compton-thick absorption, 
with a hydrogen column density ^ 10^^ cm~^. After corrected the absorption, the 2-10 keV 
luminosity becomes 10^^^^-^ erg s^^ at a distance of 16 Mpc. It is 2-3 orders of magnitude 
higher than that observed in the MEGS band (below 10 kev). We also detected a strong Fe-K 
line at 6.4^02 keV with an equivalent width of 2.4j^^ g keV, which is consistent with the heavy 
absorption. 
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NGG 3079 is an edge-on spiral galaxy at our 
vicinity (16 Mpc at Hq = 75 km s^^ Mpc^^). It 
exhibits AGN and/or starburst activities in vari- 
ous wave band. It is known as an outflow galaxy in 
radio (Duric & Seaquist 1988), (Veilleux et al. 
1994) and X-ray (Pietsch et al. 1998). In optical, 
it is classified as a LINER (low-ionization nuclear 
emission region) or a Seyfert 2, but whether the 
He line is broad or not is in controversy (Ho et al. 
1997). In far infrared, it is luminous and shows 
an extended emission, indicating a dominant con- 
tribution of the starburst activity (Perez Garcia 
et al. 2000). In radio, it is a bright H2 0-maser 
source (Nakai et al. 1995; Trotter et al. 1998; 
Sawada-Satoh et al. 2000). 

In X-ray, using the ROSAT and ASCA obser- 
vatories, Dahlem et al. (1998) and Ptak et al. 



^Institute of Space and Astronautical Science, Yoshin- 
odai, Sagamihara, Kanagawa, 229-8510, Japan 

■^Department of Physical Science, Graduate School 
of Science, Hiroshima University, 1-3-1 Kagamiyama, 
Higashi-Hiroshima, Hiroshima, 739-8526, Japan 

^Nobeyama Radio Observatory, Minamimaki-mura, 
Minamisaku-gun, Nagano, 384-1305, Japan 

^Department of Astronomy, School of Science, Univer- 
sity of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-0033, 
Japan 



(1999) reported that the NGG 3079 emission is 
dominated by a thin-thermal component and a 
power-law component below and above 2 keV, 
respectively. The former is attributed to the hot 
gas in the host galaxy. The latter with 0.4-10 keV 
luminosity of '--^ 3 x 10^° erg s~^ can be attributed 
to an assembly of X-ray binaries. Super-nova rem- 
nants and hot gas due to the starburst activity can 
also contribute to the latter component. There- 
fore, there is little room for the AGN emission in 
the ASCA spectrum. Due to poor photon statis- 
tics, the presence of Fe-K emission line have not 
been constrained with ASCA. 

We inferred that the faintness of the NGG 3079 
nucleus in X-ray is due to complete obscuration, 
and performed a BeppoSAX observation of NGG 
3079. 

2. Observation and Data reduction 

We observed NGG 3079 with BeppoSAX on 
26-27 May 2000. The exposure time of LEGS, 
MEGS2, MEGS3 and PDS are 13.6, 44.5, 44.5 
and 21.0 ksec, respectively. The HPGSPG was 
not operated. Figure 1 shows the LEGS and 
MEGS2-h3 images of NGG 3079. The brightness 
peak coincides in position with the radio/optical 
nucleus. We utihzed the LEGS and MEGS2-h3 
spectra supplied by the BeppoSAX Science Data 
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Center (SDC), which arc integrated 2' around 
the brightness peak. We utiUzed the blank-sky 
files supplied by SDC to subtract the LECS and 
MECS2+3 backgrounds. Wc also utilized the PDS 
spectrum supplied by SDC. We rebinned the spec- 
tra to have grater than 15 counts per bin. Net 
count rate of the LECS, MECS2-h3 and PDS 
are (2.9 ± 0.5) x lO'^, (4.5 ± 0.4) x lO'^ and 
(1.8 ± 0.3) X 10~^ count s~^, respectively. 

3. Spectra 

Figure 2 shows the LECS, MECS2-h3 and PDS 
spectra of NGC 3079. It is clear that the data 
requires at least three components: a featureless 
continuum in the LECS and MECS band, a heav- 
ily absorbed continuum in the PDS band, and a 
strong emission line in the MECS band. On the 
other hand, the soft thin-thermal components that 
detected in the ROS AT and ASCA spectra are not 
statistically significant in the BeppoSAX spectra, 
due to poor statistics at low energy. In this letter, 
we are not interested in the soft components, so 
that we deal here with the MECS and PDS spec- 
tra only above 2 keV, where the soft component 
is negligible. Also, we will not examine intrinsic 
absorption for the MECS-band continuum. 

3.1. The MECS spectrum 

The MECS continuum above 2 keV is well rep- 
resented with a power-law having a photon in- 
dex of ~ 2. The 2-10 keV luminosity becomes 
~ 1 X 10^*^ erg s~^. These are consistent with the 
ASCA results within the model micertainty. The 
structure around 6 keV is well reproduced with 
a Gaussian. Figure 3 shows confidence contours 
of the Gaussian component. The line; center en- 
ergy of 6.41q;| keV is in good agreement with the 
K-fluorescent line from cold iron. As shown in 
Model (a) in Table 1, the line width is not re- 
solved with the MECS energy resolution. There- 
fore, hereafter, we assume that the Fe-K line is 
narrow (Model b). The equivalent width against 
the power-law component becomes 2A~!il'^ keV in 
the 90% confidence level. It is consistent with 
the ASCA loose upper limit, although no evidence 
of Fe-K line was observed in the ASCA spectra 
(Dahlem et al. 1998; Ptak et al. 1999). Such a 
large equivalent width cannot be observed if the 
continuum emission from the nucleus is directly 



seen (e.g. Makishima et al. 1986). Consequently, 
the strong emission line suggests that the direct 
emission is completely blocked. 

3.2. The PDS spectrum 

Indeed, the PDS spectrum exceeds two orders 
of magnitude over the extrapolation of the best- 
fit Model of the MECS spectrum. Althogh the 
PDS has a large beam size (1.4° FWHM), we can 
reasonably conclude that the contribution of con- 
tamination sources to the excess is small, as fol- 
lows. In the MECS image (Figure lb; 1° circular 
FOV), there are no sources bright enough to ex- 
plain the PDS excess. For example, the observed 
fiux of the second brightest source in the MECS 
image, QSO 0957-5608, is ~ 40% of that of NGC 
3079, in the MECS band. Similarly, the sources 
in the ROSAT/PSPC image (2° circular FOV) are 
relatively faint. The 90% upper limits to the 20- 
100 keV flux of QSO 0957-5608 in the MECS and 
the other sources in the PSPC are estimated to be 
2 X 10~^^ and 3 x 10~^^ erg cm~^ s~-^, respectively. 
These arc two or one order of magnitudes fainter 
than the PDS flux in the 20-100 keV band. Ac- 
cordingly, the PDS excess is attributed to a heavily 
absorbed emission of the NGC 3079 nucleus, un- 
less there are serendipitous Compton-thick sources 
in the PDS field-of-view. 

Therefore, we fitted the PDS spectrum with 
an absorbed power-law model. Figure 4 shows 
confidence contours of the power-law component. 
The photon index is consistent with those seen in 
Seyfert galaxies, although it is not strictly con- 
strained (Model c). Therefore, hereafter, we fix 
the photon index at 1.9, a canonical value of 
Seyfert galaxies (Model d). With the fixed photon 
index, a heavy absorption with a column density 
larger than 3 x 10^** cm~^ is required. After cor- 
rected the absorption, the 2-10 keV and 20-100 
keV luminosities become ~ 1.2 and ~ 1.5 x 10^^ 
erg s~^, respectively. These are two orders of mag- 
nitude higher than that observed in the MECS 
band. 

In the spectral fitting, we considered only the 
photoelectric absorption, using the "wabs" model 
in XSPEC. In such a Compton-thick environment, 
however, the Compton scattering becomes as im- 
portant as the photoelectric absorption. As inves- 
tigated in Matt et al. (1999), the "wabs" model 
underestimates the flux 0.5-1 order of magnitude 
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in the column density range of lO^'*'^ cm ^. 
Considering the absorption due to the Compton 
scattering, the 2-10 keV and 20-100 keV lumi- 
nosities of NGC 3079 are estimated to be lO^'^-^s 
erg s~^. 

3.3. Combined fit of the MECS and PDS 

spectra 

Using the three-component model, we per- 
formed combined fits of the MECS and PDS spec- 
tra. As shown in Table 1 , the results arc similar to 
those of the individual fits, except for an improve- 
ment of the lower limit to the absorption column 
density. Model (f) in Table 1 shows our best re- 
sult. We superposed the best-fit model with each 
component on the spectra in Figure 2. The resid- 
ual on the LEGS band is roughly consistent with 
the soft component measured with ROSAT and 
ASCA. 

Starting from the best-fit model, we examined 
additional structures in the NGC 3079 spectrum. 
First, we added the "pexrav" model (Magdziarz 
& Zdziarski 1995) to the best-fit model, to study 
the Gompton-reflection hump. We fixed the pho- 
ton index of the reflection component at 1.9. The 
data prefers no reflection with an upper limit to 
the reflection efficiency of 23%. It is similar to the 
situation seen in the other Compton-thick AGNs, 
such as M51 (Fukazawa et al. 2001) and NGC 
4945 (Guainazzi et al. 2000). The smaU reflection 
efficiency indicates that not only the reflection pro- 
cess but also the transmission process contributes 
to the strong Fe-K fine of NGC 3079. Secondly, we 
applied a high-energy exponential cutoff to the ab- 
sorbed power- law component of the best-fit model. 
However, the cutoff energy was not restricted well 
with a loose lower limit of 16 keV. 

4. Variability 

In Figure 5, we show the LEGS, MECS2-h3 and 
PDS short-term light curves of NGC 3079. In the 
LEGS and MECS band, no short-term variabil- 
ity was detected during one-day observation. It 
is consistent with the spectral result that the di- 
rect component of the NGC 3079 nucleus is not 
observed below 10 keV. In the PDS band, it is dif- 
ficult to constrain the short-term variability, due 
to the poor photon statistic. 

To study long-term variability, we fitted the 



ASCA spectra above 2 keV in the same way as 
Model (b). We fixed the photon index of the 
power-law component at the best-fit value of the 
BeppoSAX data. We are not concerned with the 
absorption. The 2-10 keV flux becomes (5 ± 2) x 
10~^^ erg cm~^ s~^, so that no long-term varia- 
tion was detected during seven years. It is in con- 
trast to the three-times long-term variation of the 
M51 nucleus (Fukazawa et al. 2001). Therefore, it 
might suggest that the NGC 3079 emission below 
10 keV is mainly due to the host galaxy emission 
rather than the scattered and/or reflected compo- 
nent of the nuclear emission. 

5. Discussion 

The main result of the BeppoSAX observation 
of NGC 3079 is the detection of a heavily absorbed 
continuum and a strong Fe-K emission line. 

The column density of the absorbed continuum 
and the equivalent width of the Fe-K line allow us 
to estimate the geometry of the system. Adopting 
the geometry in Ghisellini et al. (1994), the ob- 
tained values of the column density and equivalent 
width indicate a high inclination angle of 60-80 de- 
gree, corresponding to the dominant contribution 
of the transmission component of the Fe-K emis- 
sion line. It is consistent with the no detection of 
the Compton hump in the BeppoSAX spectra. 

Using the absorption-corrected 2-10 keV lumi- 
nosity, we can estimate how active the NGC 3079 
nucleus is. Trotter et al. (1998) estimated that the 
binding mass of the inner parsec of NGC 3079 is 
^ 10^ Mq by the H2O maser kinematics. Accord- 
ingly, the X-ray emission of NGC 3079 is 1-10% 
of the Eddington luminosity. It is as active as the 
Seyfert galaxies and quasars. 

We can also examine whether the AGN or star- 
burst activity dominates the NGC 3079 energy 
output by comparing the 2-10 keV luminosity 
with the infrared luminosity. The infrared lumi- 
nosity of NGC 3079 becomes 1.5 x 10^^ erg s'^ 
(the NED data base). Assuming that the 2-10 
keV luminosity is ~ 3% of the bolometric lumi- 
nosity, as typical for quasars (Elvis et al. 1994), 
the bolometric luminosity associated to the AGN 
is estimated to be as luminous as the observed in- 
frared luminosity, hence the AGN and starburst 
activities contribute to the energy output of NGC 
3079 in the same order of magnitude. 
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Fig. 1. — (a) The LEGS image supplied by 
the BeppoSAX Science Data Center. Shown af- 
ter Gaussian smoothing with 0'.5 sigma. Con- 
tours correspond to 0.04, 0.08, 0.12 and 0.16 
counts/pixel. (b) Same as (a), but for the 
MECS2+3. Contours correspond to 0.2, 0.3, 0.4, 
0.5, 0.6 and 0.7 counts/pixel. The brightest peak 
coincide in position with NGC 3079. The second 
brightest peak 12' north- northwest of NGC 3079 
is a gravitationally-lensed quasar, QSO 0957-5608. 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v5.0. 
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Fig. 2.— (a) Crosses show the LEGS, MECS2+3 
and PDS spectra of NGC 3079. The best-fit model 
above 2 keV (Model f ) with each model component 
is superposed as histograms, (b) Residual of the 
best-fit model. Residuals below 2 keV may be at- 
tributed to the soft components observed in the 
ROSAT and ASCA spectra. 
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Fig. 3.— Two-parameter 68%, 90% and 99% con- 
fidence contours of the Gaussian component, (a) 
Confidence contours of the center energy and line 
width of the Fe-K line (Model a), (b) Confidence 
contours of the center energy and normalization. 
The dotted and solid contours correspond to the 
free line width (Model a) and that fixed at eV 
(Model b), respectively. 
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Fig. 4. - Two-paramctcr 68%, 90% and 99% con- 
fidence contours of the absorbed power-law com- 
ponent. Only the photoelectric absorption is con- 
sidered, (a) Confidence contours of the photon 
index and absorption column density. The dotted 
and solid contours correspond to the result of the 
PDS fit (Model c) and the combined fit (Model 
e), respectively, (b) Confidence contours of the 
photon index and absorption-corrected flux deter- 
mined by the combined fit of the MECS and PDS 
(Model e). The solid and dotted contours corre- 
spond to the 2-10 keV flux and 20-100 keV flux, 
respectively. 
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Fig. 5.— The LECS (a), MECS2-(-3 (b) and PDS 
(c) light curves of NGC 3079. 
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Table 1 

Best-fit parameters of the NGC 3097 spectra. 



Model 




(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


Detector 




MECS 


MECS 


PDS 


PDS 


MECS+PDS 


MECS+PDS 


Power- law 


Photon index 


2 04+°-*'' 


1 74+0.55 
^- '^-0.49 






^- '^-0.49 


1 7c, + 0.55 
'^-0.48 




2-10 kcV fiux^' 


o o + l.l 
•'■■'-1.0 


3.7 ± 0.8 






3.7 ± 0.8 


3.7 ± 0.8 


Power- law 


Photon index 






2 7(1 + 4.60 
'U-2.62 


1.9 (fix) 




1.9 (fix) 




2-10 keV flux^ 










4.915:; 






20-100 keV flux^ 
















Absorption*^ 






15.2l-:| 




i5.8i«:? 


lo.ol^-^ 


Gaussian 


Energy [keV] 


6 41+"-28 


"■■'''-0.16 






6.39j:°;i* 


6.39j:°;i* 




Line width [keV] 


44+"'^'5 

^■^^-0.44 


(fix) 






(fix) 


(fix) 




Normalization^ 




8.0 ± 4.1 








8.0 ±4.1 



X^'/d.o.f. 0.86/4 1.98/5 0.10/1 0.29/2 2.08/6 2.30/7 



0.215 0.396 0.098 0.145 0.347 0.328 



*In unit of 10~^^ erg cm~^. 

''In unit of 10"" erg s"' cm^^. 

'^In unit of lO"''^ erg s^^ cm^^. Shown after absorption correction. 

^Hydrogen column density A^h in unit of 10^^ cm~^. Only the photoelectric absorption is considered. 
^Photons in the line in unit of 10~^ cm~^ . 

Note. — Errors refer to single-parameter 90% confidence limits. The PDS/MECS normalization factor is fixed at 
0.86. 



7 



